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Executive  Summarv 


Utilizing  a  multipie-scale  decomposition  of  the  Navier-Stokes  equations,  it  has  been 
demonstrated  that  non-equilibrium  fine-scale  stress  modeling  is  sufficient  to  generate 
linearly  unstable  large  scale  vortical  disturbances.  By  assembling  various  of  these  modes, 
the  (short-time)  characteristics  of  coherent  structure  may  be  simulated  and,  eventually 
solved  via  an  initial  value  problem. 


Discussion 


The  present  research  program  seeks  to  make  fundamental  contributions  to  the  under¬ 
standing  of  wall-bounded  turbulent  shear  flows  in  constant  and  (spatially/temporally) 
variable  pressure  conditions.  Both  the  advection  and  (Lagrangian)  evolution  in  the  advec- 
ted  frame  of  reference  of  large  scale  coherent  structures  are  focused  upon.  As  a  conse¬ 
quence,  basic  deterministic  vortical  dynamics  results  as  well  as  modeling  implications  are 
sought.  In  connection  with  the  modeling,  the  deterministic  simulations  will  eventually  be 
ensembled  utilizing  a  pdf  representation  of  the  burst  occurrence,  along  with  a  temporal 
chaos  assumption  at  some  'upstream'  location. 

This  work  has  a  foothold  in  both  the  relatively  distant  past  and  current  state-of-the-art 
analysis/simulation  methodologies,  with  respect  to  the  historical  references,  the  research 
relies  heavily  on  the  following  developments: 

I.  Liepmann's  'turbular'  fluid  concept,  and 

II.  Landahl's  wave-guide  model  of  turbulence. 

If  a  thumbnail  synopsis  of  the  work  was  required  it  would  be  "to  investigate  the  merging 
of  I  and  II".  That  is,  to  extend  Landahl's  weakly  non-linear  wave-like  model  of  large- 
scale  turbulence  by  accounting  for  the  complex  apparent  constitution  of  the  medium 
induced  by  the  small  scale  processes.  That  such  an  approach  might  be  successful  is 
supported  by  comparisons  between  the  measured  and  predicted  Reynolds  stresses  from 
Crr-Scmm.erfeid  computations,  presented  by  Kim,  Kline  and  Reynolds. 
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Comparison  of  the  form  of  the  Reynolds  stress  measured  from  H,  bubble  data 
during  a  bursting  period  with  calculations  using  two-dimensional  steady  linearized  equa¬ 
tions  (Orr-Sommerfeld).  Scales  are  arbitrary.  (5,  Burst  period  data,  t  =  0-6-7  sec; - , 

two-dimensional  lines  theory. 


The  modern  analysis  foothold  of  the  research  involves  Large  Eddy  Simulation  (LES),  in 
which  the  filtered  equations  of  motion  are  used  to  resolve  deterministically  the  largest 
scales  of  motion  while  the  unresolved  sub-grid  scale  (SGS)  processes  are  modeled.  These 
SGS  processes  are  the  rigorous  theoretical  equivalent  of  Liepmann's  'turbular'  fluid. 
Unfortunately,  although  LES  shows  promise  of  being  managed  and  maturing  into  a  reliable 
predictive  tool,  it  is  often  an  un-controlled  combination  of  interacting  numericai/fiitering/ 
modeling  difficulties  and  cannot  yet  be  assumed  to  provide  stand-alone  'truth'  data.  In 
fact,  even  direct  num erica!  simulations  (DNS),  which  resolve  all  scales  of  motion,  are 
regularly  beset  with  confusion  regarding  the  true  solution  in  low  Reynolds  number  turbu¬ 
lent  flows. 
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For  this  reason,  the  research  sought  to  isolate  the  essence  of  LES,  while  avoiding  some 
of  the  problems  presently  associated  with  it.  The  multiple-scale  nature  of  the  flow  has 
been  confronted  explicitly  yet  idealized  by  empioying  the  discrete/disparate  scale  depen¬ 
dent  variable  decomposition  of  Mcllo-Christensen  and  (for  non-steady  flows)  the  coherent- 
incoherent  element  decomposition  popularized  by  Hussain  &  Reynolds  and  also  used  by 
J.T.C.  Liu.  Utilizing  the  merged  (Mollo-Christensen/Hussain  &  Reynolds)  decomposition 
(see  Attachment  1.),  the  interactions  and  transitions  between: 

*  Large-scale  temporally  coherent, 

*  Large-scale  temporally  incoherent,  and 

*  Small-scale  temporally  incoherent 

elements  may  be  focused  upon.*  Clearly,  in  this  formulation  the  large-scale  temporally 
incoherent  element  is  the  experimentally-observed  (spatially)  coherent  structure. 

As  a  result  of  this  approach  to  implementing  the  wave-guide  model  for  a  'turbular'  fluid, 
most  of  the  numerical  and  filtering  problems  of  LES  may  be  eliminated,  while  addressing 
the  small-scale  process  modeling,  as  it  affects  the  large-scale  vortical  dynamics.  To  guide 
this  modeling,  the  full  moment  equation  formulation  of  the  unresolved  processes  developed 
by  Deardorff  has  been  used.  In  so  doing,  current  algebraic/equilibrium/isotropic  approxi¬ 
mations  may  be  shewn  to  fail  in  certain  situations.  The  work  is,  therefore,  focusing  on 
the  need  for  non-equilibrium  and  anisotropic  unresolved  process  models,  relative  to  the 
prediction  of  large-s^ale  turbulent  structures. 

*The  utility  of  ''compounded)  signal-bacKground  decompositions  has  also  been  evaluated 
but,  due  to  various  m,atrematical  dif Acuities,  has  not  proven  as  susceptible  to  analysis. 


In  the  course  of  computing  weakly  non-iinear  large-scale  structure  with  various  small- 
scale  process  models,  certain  specific  deterministic  vortical  events  are  being  examined 
for  their  relevance  to  turbulence.  With  regard  to  the  physics,  the  burst  propagation  and 
evolution  of  the  structures  are  the  focus  of  interest.  At  that  point  when  results  are 
applied  to  modeling,  the  statistical  occurrence  shall  be  confronted.  Interestingly,  however, 
the  statistical  occurrence  probability  distribution  function  equation  is  generically  identical 
to  the  (L)  propagation  equation  being  studied:* 

t  c  ’  x  c  f 

where: 


Uc  =  U/H 

h  =  shape  factor,  and 
P  =  generalized  pressure  gradient. 

It  has  been  found  that,  utilizing  uc  ,  various  experimentally-derived  propagation  speed 
characterizations  may  be  understood  theoretically.  Results  for  a  given  Uc  ,  when  'un¬ 
wound'  utilizing  the  resultant  H  evolution,  imply  the  velocity  field  which  sustains  the 
flow.  In  addition,  the  dynamics  of  the  superlayer  and  interface  between  the  rotational 
and  irrotational  flow  is  dependent  uoon: 


*In  addition,  this  is  the  form  of  the  multiple-scale  perturbation  amplitude  evolution 
equation  associated  with  the:  1)  presence  of  a  coherent  non-steady  component,  2)  effect 
of  residual  nor-unearities,  and  3)  effect  of  large  scale  incoherent  on  small  scale  incohe¬ 
rent.  structures.  These  effects  are  essential  to  the  understanding  of  turbulence  perpetua- 


with  entrainment  also  proportional  to  the  amplitude  of  the  interfacial  oscillations.  The 
interaction/synchronization  between  this  outer  flew  phenomenon  and  the  other  flowfield 
instabilities  to  be  computed  is  a  primary  focus  of  the  work.  The  pressure  field  is  clearly 
a  key  diagnostic  variable  for  this  phenomenon. 

Throughout  this  work,  attention  will  be  focused  on  DL/Dt,  the  flow  normal  to  the  wall. 
Traditionally  in  laminar  boundary  layer  analysis,  the  vertical  velocity  'goes  along  for 
the  ride',  in  turbulent  flow,  the  unsteady  vertical  velocity  as  a  source  of  curvature  may 
be  a  driver,  although  non-parallel  mean  flow  effects  have  not  yet  been  shown  to  be  crucial. 

Note  that,  defining  P_1  as  a  time  constant,  t,  this  equation  may  be  written: 

t  DL/Dt  +  L  =  , 

eq. 

where: 


D/Dt  e  3/St  +  U  S/Sx,  and 
c  ’ 


L  =  TUCf  . 
eq .  c  f 


Note  that,  in  this  form,  the  equation  expresses  a  first  order  lag  between  L  and  its 
equilibrium  value,  in  advected  coordinates.  A  representative  solution  is  presented  here, 
which  illustrates  various  staoes  of  exponential  growth,  both  spatial  a^d  temporal. 


TYPICAL  l  evolution  computation 


with  regard  to  evolution  of  the  turbulent  burst  structures,  heavy  reliance  is  placed  on 
the  flcwfield  visualizations  of  Blackwelder;  Kline;  Wiilmarth,  etc.  Inasmuch  as  these 
results  gnerally  involve  streakline  patterns,  the  computational  work  is  presented  in  that 
form,  also.  Unfortunately,  the  pathline/streakiine  visualization  suffers  from  two  flaws: 


1.  Given  a  streakline,  a  unique  inference  of  the 
velocity  field  is  usually  not  possible,  and 

2.  Streaklines  can  assume  a  complex  appearance 
in  very  simple  flows  (tne  work  of  Saric  is 
illuminating  here). 


In  spite  of  these  cautions,  tne  comparison  between  computations  and  experiment  in  this 
format  is  a  crucial  test  of  the  research  hypotheses.  A  sample  pathline/streakiine  plot  is 
presented  here,  however,  extensive  use  of  this  tool  awaits  the  computation  of  the  dom¬ 
inant  unstable  modes  of  oscillation  in  the  shear  flow. 


; 

I 
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Fundamentally,  the  view  of  turbulence  adopted  here  is  relatively  traditional,  following 
along  the  lines  of  Willmarth.  Specifically,  we  assume  that  the  low  wavenumber  pressure 
field  measured  at  the  wall  results  from  the  large-scale  motions  in  the  outer  part  of  the 
flow,  and  that  these  are  triggered  by  the  bursting  in  the  buffer  layer,  which  is  caused 
in  turn  oy  secondary  instabilities  growing  on  inflectionary  profiles  induced  by  longitudinal 
vortices  in  the  sub-  and  buffer  layers. 

In  the  weakly  non-linear  formulation  which  we  are  pursuing  (prior  to  utilizing  DNS/LES 
tools)  several  of  these  physical  observations  are  (or  may  be)  related  to  classical  mathe¬ 
matical  features  of  the  Fourier/Laplace  transformed  normal  mode  vorticity  equations 
w'hich  are  being  solved.  (This  dispersion  relation/vertical  structure  distribution  calculation 
exercise  is  the  first  step  in  a  general  initial  value  problem  computation.) 
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FOURIER  AND  LAPLACE  TRANSFORMED  3-D  VORTICITY  EQUATIONS 
LINEARIZED  LARGE  SCALE  MOTIONS 


|(u-c)  +  i  [(Re-1+ 

j(u-c)  +  i  [(Re_1+ 
j(u-c)  ♦  i  [{Re-1+ 


e) (D2-k2)  +  (D2e)  -{(Du)  +1  {De)(2D2-k2)|u3 

-  |  {(De)D  -  (2D2e)}u2 
c) ( D 2 -k 2 )  +  (Dt)D]|w2  •  -  |  (Du)^ 
e)  (D2-k2)  +  (D2e)  ]|w3  •  {(Du)  +1  (De) (2D2-k2)} ^ 

+  {(De)D  -  (2D2e)}u2 


+  {-  ^  [D2u  +  (Du)D]|u2 
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First,  the  inflectional  profile  should  be  computable  from  the  U'.  evolution  equa¬ 
tions,  given  the  appropriate  r'.j  model,  which  reflects  sub-layer  structures.  This  brings 
into  question  the  role  of  critical  level  processes,  as  in  laminar  flow,  and  their  relationship 
to  the  observed  free  shear  layer  instability.  In  addition,  the  (Benney-Gustavson)  vertical 
velocity-vertical  vorticity  resonance  (or  near  resonance)  may  be  associated  with  the 
generation  of  large  vertical  structures  which  deflect  laterally  the  vertically-sheared  flow, 
thus  generating  streamwise  vorticity*:  _ 


(u-c)  4  A  [(Re*1  +  e )  (D2  -k 2  )  +  (De)D]>  =  -  -  (Du)u„ 

a  I  l  a  2. 


The  Klebanoff  hairpin  vortex  and  the  horseshoe  vortex  arising  in  local  3-D  unsteady 
separation/re-attachment  may  be  related  to  this  phenomenon. 


30  SINGULAR  POINT 


N,  ,  NODAL  POINT  OF  ATTACHMENT 


v'VVrGf5 


v  .N, ,  NODAL  POINT  OF  SEPARATION 


* 


VERTICAL 

PLANE 


D-vidi'g  surface*  formed  from  combination*  of  tu,  nod*l  point  of  attachment  and 
saddle  point,  -b)  nuda!  point  of  separation  and  saddle  point 
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Although  not  in  vogue  for  many  years,  the  effect  of  local  transient  adverse  pressure 
gradients  and  unsteady  separations  was  espoused  by  Prandtl  to  be  an  essential  aspect 
of  the  nature  of  turbulence.  The  generation  of  these  vortical  structures  and  the  presence 
of  favorable  and  adverse  pressure  fields  also  clearly  indicates  the  possibility  for  encoun¬ 
tering  classical  vortex  breakdown  phenomena.  In  its  discussion,  characterizing  the  helicity 
of  the  flow  is  also  becoming  a  useful  tool. 

Finaily,  work  has  exposed  the  role  of  spatiaily-variable  apparent  viscosities  on  the 
generation  of  streamwise  vorticity.  This  effect  is  reminiscent  of  the  formation  of  Lang¬ 
muir  vortices  and  was  carried  out  to  correct  an  error  in  the  literature,  as  well  as  to 
expose  the  effect  of  strong  vertical  variations  of  viscosity  on  the  generation  of  3-D 
vcrticity. 

In  order  to  predict  such  vortical  events,  or  the  tendency  toward  them,  the  weakly  non¬ 
linear  equations  of  motion  are  being  solved  with  various  un-resolved  process  models. 

The  fundamental  problems  with  this  aspect  of  the  work  involve  the  non-uniqueness  of 
such  modeling  and  the  lack  of  data  on  small-scale  turbulence,  conditionally-sampled  to 
reflect  its  dependence  on  the  burst  events.  Therefore,  as  systematically  as  possible,  the 
work  is  tediously  surveying  the  effect  of  non-equilibrium,  anisotropy  and  non-linearity 
in  the  un-resolved  process  model  on  the  large-scale  structure  prediction.  At  this  point, 
due  to  the  importance  of  the  effect  on  turbulence  modeling  in  general  and,  also,  for 
visco-elastic  m.edia*  (see  Crow;  Davis;  Betchov  and  Criminals),  the  effect  of  non-equili- 
crium  is  being  concentrated  upon.  Based  on  a  recurrent  Symmetry  of  various  turbulent 
process  modeling,  the  first-order  lag  model  has  also  been  adopted  here: 

•The  small-scale  turbulence  field  is  known  to  behave  in  such  fashion  also. 
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Various  candidate  non-linear  and  anisotropic  'stress'  models  are  also  being  analyzed. 
Unfortunately,  the  non-uniqueness  problem  is  even  more  prevalent  here.  The  work  is, 
therefore,  attempting  to  determine  whether  or  not  small  scale  structure  non-equilibrium 
is  an  essential  aspect  of  the  turbulent  burst  dynamics  processes. 


The  analysis  is  able  to  handle  arbitrary  velocity  profiles  as  a  function  of  space  and  time. 
However,  relative  to  TBL  structures,  the  data  to  validate  a  large-eddy  model  is  concen¬ 
trated  in  simpler  flows.  Therefore,  until  the  essential  features  of  the  un-resolved  process 
modeling  are  clarified,  classical  data  shall  be  the  focus  of  simulation.  In  fact,  to  supple¬ 
ment  the  constant  pressure  TBL  simulations  which  have  been  on-going,  the  transitional 
flow  of  Schubauer  and  Klebanoff  is  being  addressed  to  prove  the  model  in  a  situation  in 
which  a  (streamwise)  change  in  instability  processes  is  occurring,  as  the  mean  vorticity 
distribution  changes  and  unresolved  process  effects  arise. 


Typical  results  for  an  unstable  mode  of  the  predicted  disturbance  vorticity  are  presen- 
tea  here  for  a  (first-order  lag)  unresolved  process  model.  Various  parametric  sweeps 
through  the  burst: 


1.  Frequency, 

2.  Span-wise  wave  number, 

3.  Complex  stream-wise  wave  number, 

4.  Nor, -equilibrium  time-cc-nst art,  ar.g 

5.  Unresolved  process  model  constants. 


are  being  carried  out  for  these  transitionai  profiles  to  isolate  the  essential  ingredients 
for  a  small-scale  turbulence  model  which  allows  the  observed  large  scale  structure 
streaklines  and  pathlines  to  be  simulated. 


I 

I 
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The  conditions  for  which  this  outer  flow  disturbance  was  found  are: 


8  =  0.1 
to  =  2.0 
e  =  100 

Te  =  1.0 

Re  =  4.4369  x  10" 


(Span-wise  wavenumber) 

(Frequency) 

(Peak  small  scale  stress  divided  by  v) 
(Non-equilibrium  time-constant) 


where  all  variables  are  non-dimensionaiized  by  6*  and  U^.  The  streamwise  wavenumber 


solution  is: 


a  =  .162  -  0.0373. 


Note  that  the  angle  of  the  disturbance  in  the  horizontal  plane  conforms  with  experimental 
observations  also.  For  this  case  the  amplitude  distribution  of  the  small-scale  stresses  was 
fit  to  the  function  y2  e~Ay  in  order  to  simulate  traditional  experimental  data  such  as 
those  in  the  attachment  here. 


Equally  important,  however,  is  the  population/character  of  the  solutions.  A  typical  plot 
of  the  boundary  condition  error  (wall  normal  velocity)  is  shown  here.  Note  that  the  various 
error  spikes  always  neighbor  zero  error  points,  i.e.  solutions.  ror  the  disturbance  solution 
presented,  the  coarse  resolution  and  fine  resolution  error  fields  are  also  presented. 
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Conclusions 


I 


During  the  forthcoming  year  of  work,  additional  variable  pressure  velocity  fields  shall 
be  addressed  as  the  quality  of  the  unresolved  process  modeling  warrants.  In  addition, 
the  research  hypothesis  (that  the  merging  of  the  wave-guide  model  of  turbulence  with 
an  appropriate  complex  small  scale  turbulence  constitutive  model  will  clarify  some  of 
the  turbulence  vortical  dynamics  and  the  modeling  thereof)  shall  be  thoroughly  tested, 
if  the  weakly  non-linear  modeling  is  successful,  the  SGS  modeling  implications  for  LES 
will  be  significant.  If  the  wave-guide  model  is  not  susceptible  to  "improvement",  it  will 
be  clear  that  full  non-linearity  is  required  to  simulate  even  the  initial  stages  of  the  burst 
process  and  that  the  un-resolved  process  modeling  is  not  critical,  as  long  as  it  reflects 
reasonably  the  scales  of  motion,  in  such  an  event,  the  wave-guide  model  would  be  viewed 
as  a  tool  to  understand  the  nature  of  an  existing  fluctuation  field  rather  than  its  incep¬ 
tion  and  evolution. 
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A  FOUR-ELEMENT  DECOMPOSITION  OF  THE  NAVIER-STOKES  EQUATIONS 
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Decompose  the  flow  Into  four  components  by: 
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Incoherent  Large  Scale  Flow: 
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Incoherent  Small  Scale  Flow 
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M?an-ve loclty  profiles  through  transition  region. 

80  feet  per  second*  free-stream  turbulence,  0.05  percent 


'iLes  of  u  /U^  through  transition  region.  *  80  feet 
second.’,  free-streeun  turbulence,  0.03  percent. 


